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ABSTRACT 



This thesis presents a mixed integer program (MIP) to 
schedule sequences of classes attended by Marines at the 
Marine Corps Communication- Electronics School in order to 
qualify them for communications and electronics military 
occupational specialties. The schedule determines the 
starting dates and the number of students to enroll in each 
instance or "class" of various course types. The courses 
follow a specific sequence and many classes of a course may 
be scheduled within a fiscal year. Students attend one or 
more of the courses and may wait some time for a class of a 
subsequent course to convene. The objective of the MIP is to 
reduce the amount of delay students incur while waiting for 
classes of additional courses to start in the sequence. Due 
to the size and complexity of the model, the MIP initially 
schedules classes with a weekly resolution but then adjusts 
the starting dates to produce a daily schedule. For 1993 
data, the MIP is solved in less than 10 minutes on a desktop 
computer (80486 processor at 66Mhz and at least 64M RAM) and 
produces a schedule which has 62% less delay than the actual 
schedule for that year. 
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THESIS DISCLAIMER 

The reader is cautioned that computer programs developed 
in this research may not have been exercised for all cases of 
interest. While effort has been made, within the time 

available, to ensure that the programs are free of 

computational and logic errors, they cannot be considered 
validated. Any application of these programs without 
additional verification is at the risk of the user. 
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EXECUTIVE SUMMARY 



This thesis presents a method to improve the efficiency of 
the training of enlisted Marines in communications and 
electronics Military Occupational Specialties (MOSs) . A mixed 
integer program (MIP) is developed which schedules the courses 
provided by the Marine Corps Communication- Electronics Schools 
(MCCES) . The objective is to minimize the delay students 
incur while waiting for follow-on courses. 

To meet the annual MCCES training requirements (e.g., 
90,804 man-days of class work in 1993) , more than one instance 
of each course is usually required. Each instance of a course 
is designated a "class" and a completed schedule requires 
specifying the starting dates and the number of students to be 
enrolled in each class. The courses follow a specific 
sequence with students proceeding to subsequent courses after 
successful completion of the initial course. The number of 
students who can be scheduled for a class of a course must 
remain between a maximum and minimum for each course type. 
The number of concurrent classes of a given course must remain 
below a maximum for that course and a minimum number of days 
between starting dates of classes must be maintained. 

These constraints may cause students to incur some delay 
while waiting for classes of a subsequent course to begin. 
The goal of this thesis is to develop a schedule which 
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minimizes the student delay so that students may proceed more 
quickly to their subsequent assignments. 

The size of the MIP created by a set of 17 courses and a 
planning horizon of 236 training days is large and a solution 
cannot be obtained. When the planning horizon is scaled to a 
weekly resolution, the size of the MIP is reduced by nearly 
80% and solution times become manageable. A rescaling of the 
planning horizon back to a daily resolution along with a 
partial correction for delay caused by scaling is then 
performed. Fractional attrition of students does create 
difficulty in obtaining schedules having integer class sizes, 
but this is largely overcome by assuming classes have known 
integer values for attrition. 

Regardless of the method used to schedule classes, 
limitations in resources will impose some delay of students in 
completing training. However, an existing 1993 schedule of 
classes for a set of courses from MCCES exhibits 7,819 man- 
days of student delay, which adds 8% to the total man-days 
required for training, while a schedule developed by the MIP 
using the same 1993 data adds 2,964 man-days of delay which is 
only 3% of the total. Schedules created by the MIP can be 
obtained on a desktop computer (80486 microprocessor at 66MHz 
with at least 64M RAM) in approximately 30 minutes. 
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I . INTRODUCTION 



Every entrant into the armed services must undergo initial 
training to qualify for service in a specific military 
occupational specialty (MOS) . In the Marine Corps, the 
entrant must pass a series of one or more courses to acquire 
an MOS. This thesis presents a mixed integer programming 
(MIP) model to aid scheduling of courses and students in 
communications and electronics MOSs. The goal of the model is 
to reduce the waiting time of the students between subsequent 
courses so that training is accomplished efficiently and the 
students move to their subsequent assignments more quickly, 

A. BACKGROUND 

The Marine Corps Communication Electronic School (MCCES) 
located at the Marine Corps Air Ground Combat Training Center, 
Twentynine Palms, California, is responsible for the training 
of entry level Marines in electronic fundamentals, operational 
communications, air control/antiair warfare operations and 
communications -electronics maintenance. MCCES has developed 
specific courses of instruction designed to meet the 
requirement for MOS training in these areas. 

Headquarters Marine Corps identifies the annual training 
requirement for each course in terms of total number of 
students who must start training and then directs the school 
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to develop a schedule for the courses it provides. The 
schedule must identify the number of students in and the 
starting date of each occurrence of a course. The number of 
possible days which may be considered for starting courses are 
referred to as "training days" and excludes weekends, federal 
holidays and miscellaneous other days designated by the 
commander of MCCES. 

Each occurrence of a course is called a "class" and is 
designated by course name, a number identifying which 
occurrence it is in the fiscal year and by the fiscal year in 
which the class will be conducted. An example of the class 
designation for the first basic electronics course (EEC) 
conducted in fiscal year 1993 would be, EEC 1-93. Typically, 
more than one class for each course is necessary to satisfy 
the year's training requirements because the requirements can 
be large and the class sizes are limited. The enlisted 
assignment branch at Headquarters Marine Corps uses the 
schedule produced by MCCES when they identify and order 
students to receive MOS training after successful completion 
of other initial training. Eoth officer and enlisted 
personnel undergo training at MCCES. However, since the 
number of officers is only a small percentage of the total 
number of personnel trained at MCCES, this thesis considers 
the scheduling of classes for enlisted personnel only. 

Much of the MOS training provided by MCCES requires the 
progression of students through a sequence of courses, one at 
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a time, beginning with an elementary course followed by a 
number of specialized courses. Courses in a sequence have a 
specific ordering which must be followed, establishing a 
predecessor- successor relationship. For each course (except 
the first course) in a sequence, a predecessor course is 
uniquely defined. Therefore, all the courses a student has 
attended upon completion of training at MCCES may be 
identified by the last course taken and the unique predecessor 
property. 

Upon completion of a course, a student will be associated 
with one of three exclusive categories: 1. The student has 
completed training at MCCES. 2. The student starts a 
successor course on the next training day after completion of 
a predecessor course. 3. The student is placed in 
"inventory" awaiting additional training. 

Progression of students through a sequence of courses, one 
at a time, is a common method of training for many military 
MOSs, especially those MOSs which are highly technical. Pilot 
training, electronics maintenance and avionics are examples of 
the types of MOSs which require completion of sequential 
courses. Much of the research in this thesis is also 
applicable to such MOSs. 

B. RELATED STUDIES 

There has been much written on the subject of scheduling. 
The problem presented in this thesis has both similarities and 



3 



differences with many types of scheduling problems and a 
discussion of scheduling problems most closely related to the 
MCCES scheduling problem follows. 

MOS training of Marines is somewhat analogous to a 
frequently seen scheduling problem often referred to as 
"timetabling". The timetabling problem arises when students 
and a teacher (or teachers) must meet over a finite period of 
time, like a semester, and require resources, e.g., rooms, 
teaching aids, etc. [Ref. 1] . The nature of this 
problem is that a fixed number of classes of a course type 
must be assigned to specified blocks of time within the 
planning period. Timetabling typically does not involve the 
scheduling of courses which have a predecessor- successor 
relationship and differs from the MCCES scheduling problem in 
this regard. Additionally, a student at MCCES attends only 
one course at a time and there is no need to deconflict two or 
more courses competing for assignment to a specific time 
block. 

Attendance at one or more courses which have a specified 
sequence is similar to "job-shop" scheduling. As defined by 
Gere [Ref. 2], job-shop scheduling specifies a number 
of jobs, each comprising one or more operations to be 
performed in a specified sequence on specified machines and 
requiring certain amounts of time, to be scheduled such that 
due dates associated with each job will be met, or, failing 
this, some measure, such as the amount of time exceeding the 
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due date, is minimized. Davis [Ref. 3] lends insight 
to this broad definition by observing that, typically, the 
jobs are assumed to be one-of-a-kind orders, with simple 
(mostly serial) precedence orderings among operations, and 
requirements of only one machine for each operation. Also, it 
is generally assumed that only one machine is available for 
each operation. The Marine Corps MOS training problem differs 
from the job- shop scheduling problem in that the number of 
jobs (students) to be processed is variable between certain 
limits, may vary from machine to machine (course to course) , 
and there may be more than one machine (course) of each type 
running at a given time. 

Another similar type of scheduling problem is that of the 
"resource -constrained project scheduling" problem. As 
presented by Fisher, [Ref. 4], the resource- 
constrained project scheduling problem requires the 
determination of a series of start times of specified tasks 
which minimize some function of the task completion times. 
The tasks are constrained by predecessor- successor 
relationships expressed by a network and expend specified 
amounts of various scarce resources. The resource constrained 
project scheduling problem addresses many of the desired 
objectives for project completion including minimizing project 
duration, minimizing project cost given performance payments 
and penalties, and minimizing the consumption of critical 
resources [Ref. 5]. This type of scheduling problem 
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is related to the MCCES scheduling problem in that resources 
(instructors) are expended in accomplishing the tasks 
(courses) and the tasks have a specified predecessor- successor 
relationship. However, like the job- shop scheduling problem 
this problem does not address the variability in job (class) 
sizes from task to task (course to course) . 

C. PURPOSE AND SCOPE 

The purpose of this thesis is to reduce the amount of 
delay encountered by students waiting for additional courses. 
The thesis will identify a schedule of classes which specifies 
the starting dates for classes of each course, along with the 
number of students which must attend each class. 

The model presented in this thesis may be used for the 
scheduling of any courses provided by MCCES, but the primary 
focus is on courses which provide entry level training of 
students in communications and electronics maintenance. The 
model is not intended to replace the scheduler. The scheduler 
will still need to apply judgement to adjust the schedule to 
accommodate changes in requirements, resources, or course - 
specific data which occur mid-year. 

The number of students finishing a class of course c is 
not of principal importance because a completed schedule 
requires identification of the number of students starting a 
class of course c. The difference between the number of 
students starting and finishing a class of course c. 
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attrition, will be considered when determining an estimate of 
the number of students finishing a class of a course with 
subsequent assignment to a class of a successor course. 

D. OUTLINE 

In Chapter II, gives details of MOS training in the 
communications and electronics fields. Chapter III presents 
two mathematical formulations of a MIP for the scheduling of 
classes at MCCES. The first MIP is a basic model and the 
second is a refinement of the basic model to allow for any 
necessary violations of constraints. Chapter IV describes the 
implementation of the mathematical model and Chapter V 
presents computational results. In the final Chapter, a 
discussion of the solution to the integer program is presented 
along with recommendations for further development and 
improvement . 
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II. PROBLEM DESCRIPTION 



The problem faced here is the development of a schedule of 
classes for the courses MCCES provides. This problem involves 
establishing the number of classes of a course necessary to 
meet the published requirements ; identifying starting dates of 
each class and specifying the number of students to start each 
class . 

A. NEED FOR A SCHEDULE 

MCCES must administer 53 different courses and manage 
approximately 4000 students to meet the Marine Corps' annual 
training requirements for MOSs in the communications and 
electronics fields. A schedule to achieve this must be 
developed each year to meet changing requirements. 

The identification of the starting dates of all the 
classes of each course assists the Marines Corps assignment 
branch in matching the arrival of students to MCCES with the 
starting dates of the initial class they are to attend at 
MCCES. This policy is intended to minimize the time students 
wait for a class to start. For students having completed a 
class of a course and designated to attend a successor course, 
the schedule of classes will identify the earliest date the 
students may start the successor course. The amount of time 
the students spend waiting for a class of a successor course 
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to begin after completion of a predecessor course is a measure 
of efficiency of a schedule. A schedule of all classes of 
each course also provides a means to determine when students 
will complete training and proceed or return to duty. 

B . SCHEDULE DEVELOPMENT 

A completed schedule has many components, but the critical 
elements are the starting date of each class and number of 
students scheduled for each class. The other elements of the 
schedule, class ending date, class designation, and date 
students must report to MCCES (in the case of class of a 
course with no predecessor) , can be derived from the first two 
elements . 

The total number of students to be scheduled for training 
in each course for the next fiscal year is developed by 
Headquarters Marine Corps and is published in a Training Quota 
Memorandum (TQM) in December of each fiscal year. Class 
scheduling is currently done manually by MCCES personnel 
modifying the previous year's schedule, making any adjustments 
necessary to account for changes in the demand. This manual 
scheduling makes as few changes as possible to the existing 
schedule and thus does little to change the amount of time 
students spend waiting for a successor course to begin, or 
equivalently, change the average number of students in 
inventory. Once a schedule is completed by MCCES, it is 
forwarded to Headquarters Marine Corps for endorsement. 
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C. IMPROVING THE SCHEDULE 



At this time, there is no measure in use for the amount of 
delay created by the current scheduling method. For the 
purpose of this thesis, delay is defined as the days a student 
waits to start a successor course after completing a 
predecessor course. Each day a Marine spends waiting for a 
class of a successor course to begin is one man-day of delay. 
An improved schedule is a schedule which minimizes total man- 
days of delay over the planning horizon. Improving the 
schedule benefits the Marine Corps by enabling Marines to 
return to duty at the earliest possible time after having 
completed all required training. 

D. STRUCTURE OF THE COURSES 

Due to the related material presented in many of the 
courses at MCCES, the courses are divided into mutually 
exclusive groups. Each grouping has a unique initial course 
and all courses within a group have a particular ordering 
establishing one or more sequences of courses. Many of the 
course groups consist of a single sequence of courses having 
one or more courses in the sequence. Development of a 
schedule to reduce waste in this type of group may be trivial 
in that a class of a successor course may start as soon as 
sufficient number of students have completed classes of the 
predecessor course. 
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A far more complicated grouping of courses is one which 
has as many as 17 courses, has multiple course sequences, with 
each sequence containing up to four courses. A network 
representation of this type of course group is depicted in 
Figure 1. The MIP presented in Chapter III will handle course 
groups of this level of complexity. 





Figure 1. Directed network diagram of a course grouping 
beginning on the left with the unique initial course. Nodes 
i represent courses, arcs (i,j) indicate the course j is 
uniquely preceded by course i. Each echelon, left- to-right, 
represents a transition to successors courses from the 
prerequisite course. 
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III. MODEL DEVELOPMENT 



This Chapter develops the mathematical model for 
scheduling and explains its features. The basic model has a 
generalized network flow component identifying the flow of 
students through a sequence of courses to meet the annual 
training requirement for each of the course types. However, 
the addition of constraints such as maximum and minimum class 
size complicates the model and increases the difficulty of 
obtaining an optimal solution. 

A. BASIC MODEL 
1. Indices 

c,c" - courses 
t,T * training days 
C - set of courses in a grouping 
C" - set of courses which are successors to 
course c, C" C C 

T - { -235 , -234 , . . .235 , 236 } , set of training days 
in two consecutive fiscal years. The positive 
elements represent the training days of the 
fiscal year being planned and all other values 
represent training days in the immediately 
preceding fiscal year 
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2 . Data 



MAX, 


- upper bound on class size for course c 


MIN, 


- lower bound on class size for course c 


Lc 


- number of training days necessary to 
complete a class of course c 


TIP, 


- demand for course c in the planning period. 
The number of students required to start 
training in course c during the planning 
period 


CPTY, 


- maximum number of concurrent classes of 
course c 


INTV, 


- minimum number of days between starting 
dates for course c 


ATTR, 


- fraction of students starting course c who 
fail to complete course c 
0 s ATTR, s 1 


Yet 


- number of students in a class of course c 
beginning on day t for (1-L,) s t s 0 


Heo 


- number of students in inventory on day 0 


Zet 


- parameter equal to 1 if a class of course 
c starts on day t for (1-L,) s t s 0, and 
equal to 0 otherwise 


Ac 


- cost of a student who has completed a class 
of course c and is in inventory 
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MAXp, MINg, and L^. are obtained from the Course 
Descriptive Data (CDD) , a document which is periodically 
reviewed and updated and remains on file with the individual 
course coordinators. The value of L^. is fixed and is 
proportional to the amount and difficulty of the material 
presented in each course syllabus. The type of instruction, 
lecture or lab, also contributes to determining course length. 
MAXp is limited by size of classrooms being used, amount of 
equipment available for use in lab sessions and the number of 
instructors assigned to teach classes of course c. MIN^ 
identifies the minimum number of students in a class of course 
c which ensures resources are efficiently used. For instance, 
the scheduling of a class for one student requiring two 
instructors and a classroom wastes valuable resources. Having 
a student remain in inventory until MIN, students are 
available to start course c is usually preferred to expending 
valuable instructor (and other) resources to train a small 
number of students. 

A Training Quota Memorandum (TQM) is published once a 
year and specifies values of TIP, for each course. The values 
identified in the TQM for TIP, are expressed in terms of the 
total number of students which must be scheduled to start each 
course in the upcoming fiscal year. 

INTV, and CPTY, must be determined by persons familiar 
with day-to-day operation of each course. These values are 
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based on limitations such as classrooms being shared by other 
classes of course c, the qualifications and number of 
instructors assigned to the course, and the type of equipment 
used during practical application sessions. INTV^ is most 
closely related to the number and qualifications of 
instructors assigned to a course. Material presented in a 
course is segregated into blocks of instruction with a block 
usually having duration INTV^. Blocks of training within a 
course routinely have a specific ordering which must be 
followed. It is possible that not all instructors assigned to 
teach a course are qualified to teach all blocks of a course. 
If only one qualified instructor is available to teach the 
first block of a course, and there is a class undergoing 
instruction in the block, no additional classes of course c 
may begin until that instructor is available which will occur 
after INTV^. As another example, consider a block of training 
which requires the use of certain test equipment which is only 
available in sufficient quantity for one class. Once a class 
that is using the equipment has completed the block of 
training, another class could enter the block of training 
which uses the equipment. 

The values of are obtained from the schedule of 
classes for the fiscal year immediately preceding the fiscal 
year for which classes are being scheduled. The Y^t are only 
defined for {l-L^) s t s 0 which implies a class starting on 
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day t finishes in the fiscal year being planned or finishes on 
the day immediately prior to the first day of the fiscal year 
being planned. is the number of students who have 

completed course c and are in inventory on the day immediately 
prior to the first day of the fiscal year being planned. 
identifies a class of course c being started on day t for (1- 
Lp) s t s 0 and equals 1 when Y^^ > 0, and equals 0 otherwise. 

The cost per day of students in inventory, A^., is 
specified by the user. It may be the per day salary of the 
typical rank of the student in inventory or any other 
reasonable measure of interest. 

3. Decision Variables 

- number of students in a class of course c 
beginning on day t. y^ = 

(1-LJ s t s 0 (this handles classes which 
start in the previous fiscal year but finish in 
the current planning period or on the last day 
prior to the current planning period) 

Zpt - binary variable equaling 1 if a class of course 
c starts on day t and 0 otherwise. s z^^ 

for (1-Lj 5 t s 0 

h^t - number of students in inventory (waiting for 
the next course) on day t, having completed a 
class of course c. h^o = 
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4. Formulation 



Stated simply, the overall objective is to get 
students through a sequence of courses with as few man-days 
of delay as possible. The lengths of the individual courses 
are fixed and the demand for each of the courses is also fixed 
so the time a student spends waiting between the end of a 
course he will attend and the beginning of any subsequent 
courses he may attend must be minimized. To achieve this 
objective, a function which minimizes the total man-days of 
inventory is used in the basic MIP model. 

The mathematical description of the model is as 

follows : 

Basic Model 

minimize 

C til 



Subject to: 

Y, Yet = TIPc V c ( 1 ) 

til 

£ MAX^Z,,^ Vc,tal (2) 

y„ a Vc,tal (3) 

t 

52 ZcT ^ 1 

T=t- (INTVc-1) 
t 

z„ s CPTY^ Vc,tal (5) 

t=t- (L<--1) 
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hc(c-i)^ a Y, Vc'c * Vc,t 2 l (6) 

c' € 

0 ^ V c , t ^ 1 

0 ^ V c, t^l 

€ {0,1} V c, t ^ 1 
V c,l-L^^t^0 
h,, ^ H,, V c,l-L,^t^0 

Zct - Z„ Vc,l-L^sts0 

5. Explanation of Constraints 

Constraints (1) require that the total number of 
students scheduled to start each course over the planning 
period must meet the demand. Constraints (2) and (3) are 
variable upper bounds and variable lower bounds, respectively. 
These constraints specify that if a course is scheduled to 
begin on day t, the course must contain no fewer than MIN^ 
students and no more than MAX^ students. 

Constraints (4) stipulate that no more than one course 
of a given type may begin within a specified interval. Also, 
the number of concurrent classes of a given course type must 
remain within the capacity for that course, represented by 
constraints (5) . 

Constraints (6) require that the number of students 
scheduled to start a successor course on day t not exceed the 



18 



number of students which have completed that course's unique 
predecessor. The students scheduled for a successor course 
may be drawn from a class of a predecessor course completed on 
the day prior to the current day or may be drawn from an 
inventory of students who have completed a class of the 
predecessor course prior to the previous day. It is not 
necessary for a predecessor course to have started in the 
current planning period. 

Constraints (6) are generalized network flow balance 
constraints stated as inequalities because not all students 
are required to attend an additional course. Students only 
proceed to successor courses in accordance with a successor 
course's training requirement. The flow balance constraints 
are "generalized" because of attrition: the exact numbers of 
students who finish a class they start is not known but an 
estimate of this number is (1 - ATTR^y^f Of course, the 
number of students who finish a class is an integer but since 
ATTR^ is an estimated fraction, (1 - ATTRc)y^, will also be a 
fraction in most cases. This causes noninteger coefficients 
in constraints (6) and will likely cause noninteger class 
sizes in the schedule. The issue of noninteger class sizes is 
addressed in Chapter IV. 

Examination of the original formulation reveals the 
potential for tightening the linear programming relaxation of 
the model by adding constraints which are redundant in 
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defining the integer feasible region but remove unnecessary 
parts of the feasible region for the model's continuous 
relaxation. Let fx! denote the smallest integer greater than 
or equal to x and let denote the largest integer less than 
or equal to x. The number of classes scheduled for course c 
must be at least / MAX,"] and at most [tiPc / MIN^ J and 

these requirements can be explicitly added to the model: 

S [( TIP^/MAX^)l Vc ( 7 ) 

til 

Y, Zct ^ TIP^/MIN^)J Vc (8) 

til 

B . ELASTIC MODEL 

Elastic modeling allows for the penalized violation of a 
set of constraints and can assist in identifying constraints 
which otherwise would render a model infeasible in the classic 
sense. It is possible that a given set of data may prevent a 
solution to the model specified by the constraints. 
Constraints (1) through (8) are rigid constraints in that they 
do not allow for any violations to occur without rendering the 
model infeasible. For instance, given CPTY^ and MAX^., it may 
simply be impossible to meet the demand TIP^. 

Violation of a set of constraints may be reasonable. 
Consider an example where course c has a minimum class size of 
five. If four students are presently available to start 
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training in this course in the rigid model they must remain in 
inventory until at least one additional student is available. 
Constraints (3) may be modified to allow four students to 
start course c and prevent any further accumulation of man- 
days in inventory. It may be preferable to "pay a penalty" 
and violate the minimum class size restriction rather than 
have students remain in inventory for an extended period of 
time, especially if the system has enough capacity for an 
additional class of course c. Except when students must work 
in groups of size MIN^ in order to carry out certain 
experiments, class size restrictions are set for convenience 
and slight deviations may be acceptable. 

However, elastic versions of certain constraints of this 
model may be neither acceptable nor applicable. Constraints 
(1) , (2), (3), (4) and (5) will be elasticized to allow 
penalized violations of the demand, class size maximum, class 
size minimum, minimum starting interval of consecutive classes 
and maximum capacity of concurrent classes of course c, 
respectively. Constraints (6) must remain rigid because 
material presented in successor courses requires the mastery 
of topics presented in predecessor courses. 

Since the previous formulation of the model does not 
address the issues identified above, the following model is 
presented as a modification to the basic model. 
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1. Additional Indices 



None . 

Additional Data 

Bgj - cost of beginning a class of course c on 

day t 

MAXPENgt - penalty per student for violating maximum 
class size of course c starting on day t 
MINPEN^t - penalty per student for violating minimum 
class size of course c starting on day t 
TIPPENj - penalty per student for exceeding the 

annual training requirement of course c 
TIPPEN; - penalty per student for failing to meet 

the annual training requirement of course 
c 

FREQPENj. - penalty per class for violating the 

minimum interval for course c 
EXTRAPENg - penalty per class for violating the 

maximum number of concurrent classes of 
course c 



The data from the initial formulation, MAX^., MIN^, TIP^., 
INTVg, and CPTY, are unchanged. The parameters introduced for 
the reformulation must be specified by the user with values 
assigned by the relative importance the user places on the 
violation of the constraints to which the penalties apply and 
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relative to the inventory cost per student per day A^. These 
penalties appear as coefficients of elastic variables in an 
objective function presented later. Bg, is added to tune the 
behavior of the model and is the cost on the class scheduling 
variable 

3. Additional Variables 



Sci 


slack variable identifying the number of 
students above MIN, for a course 
starting on day t, 

0 s s„ s MAX, - MIN, 


overmax,, - 


number of students exceeding MAX, for 
course c starting on day t 


undermin,, - 


number of students under MIN, for course 
c starting on day t 


overtipc 


number of students exceeding TIP, for 
course c 


undertip. 


number of students under TIP, for course 
c 


freqvio,. 


number of extra ( >1 ) classes of course 
c starting in the interval 
[t- (INTV,-1) , t] 


Vci 


slack variable identifying the number 
( < 1 ) of classes of course c starting 
in the interval [t- (INTV,- 1) , t] (used in 
reformulation of a set of constraints) 



23 



extracl, 



number of concurrent classes of course 



c exceeding CPTY^ 

4. Elastic Formulation 

The new objective function includes the student 
inventory variable of the original objective function along 
with all of the elastic variables and their associated 
penalties, and the class scheduling variable and its new cost. 

The mathematical description of the elastic model is 
as follows: 

Elastic Model 



minimize Y, T, + E E ^ct^ct + 



C€C til C£C til 



C£C til 



5^ MINPEN^undermin^t ^ ^ MAXPEN^overmax^^ + 




Y Y FREQPEN^freqvio^t ^ YY EXTRAPEN^extracl^t 



C€C til 



ceC til 



J^TIPPENc undertip^ + TIPPENc over tip^ 



ceC 



Subject to: 



Y Vet under tip^ - overtip^ = TIP^ V c 



(El) 



Vet ^ct undermin^,^. - overmax^^ = MAX^z^,. 



Vc,t^l (E2/3) 



t 

53 ^CT = 1 V c,t = l (E4a) 

T = t- (INTVc-1) 

^ct “ ^c(t-INTV_) ^ ''^ct “ ’^c(C-l) _ ^ 

(E4b) 

- freqvio^^ + f reqvio^(^_i) =0 V c,t^2 

c 

J3 - extracl^.^ ^ CPTY^ V c,t^l (E5) 

t=t-(Lc-l) 

hc(t-D* (l-ATTRc)y<,,t-i„) i E Vc't + Vc,t2l (e6) 

C' € €*= 

Z« i f( TIP^/MAX,)1 Vc (E7) 

til 

j; U TIP^/MIN^)J Vc (E8) 

til 

0 ^ V c, t 



0 £ V c, t 

0 ^ ^ MAX^ - MIN^ V C, tkl 

z^t 6 {0, 1} V c, t 



25 



V C, l-L^^ t ^0 



y ct ^ct 

= Zct V C, 1-L^^t ^0 

5. Explanation of the Elastic Model 

The new objective function includes all elastic 
variables and their associated penalties along with the 
inventory variable and its cost and the class scheduling 
variable and its cost. A positive value of an elastic 
variable should worsen the objective function value therefore 
all elastic penalties are positive. The cost of the class 
scheduling variable is actually negative which encourages 
classes to begin. This is necessary because the form of the 
elastic constraints (E2/3) might allow y^j and overmax^t to both 
take on positive values with z^t still equal to zero, meaning 
students could be assigned to a class which is not scheduled. 

Constraints (El) are a modification of constraints (1) 
to allow a penalized deviation from the demand for course c. 
Constraints (E2/3) combine the class minimum and maximum size 
constraints into a single elastic constraint which allows 
penalized violations of the minimum and maximum values. These 
constraints are combined by inserting a bounded slack 
variable, This method of combining of linearly dependent 

constraints is described by Brown and Olson [Ref. 6] 
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and eliminates a block of constraints at the expense of adding 
an additional variable s,,,. 

Constraints (E4a) and (E4b) are a modification of 
constraints (4) to allow for a penalized violation of the 
number of classes of course c starting in the interval [t- 
(INTVg-1) , t] if the number of classes in the interval is 
greater than one. Additionally, constraints (E4a) and (E4b) 
are reformulated as “pure network” constraints for reasons to 
be discussed in Chapter IV. 

Constraints (E5) modify constraints (5) by the 
addition of the elastic variable extracl^t. The inclusion of 
extraclct allows for a penalized violation of the maximum 
permitted number of concurrent classes of course c. 

Constraints (E6) remain unchanged from the initial 
formulation. They are not elastic because the predecessor- 
successor relationship must not be violated . A breach of 
these constraints would mean students could be scheduled for 
a successor course without having completed that course's 
unique predecessor. 

Constraints (E7) and (E8) are identical to constraints 
(7) and (8) of the original model. The purpose of these 
constraints in the rigid model is to tighten the linear 
programming solution of model for a quicker solution to the 
MIP. The inclusion of these constraints in the elastic model 
accomplishes the same goal. With constraints (E7) and (E8) 
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remaining rigid, constraints (E2/3) are restricted to minor 
violations of MAX, and MIN,. 



I 

I 
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IV. MODEL IMPLEMENTATION 



The MCCES scheduling problem represented by the elastic 
model of Chapter III is solved by coding the model in GAMS 
(General Algebraic Modeling System) and solving the model 
using the X- System. GAMS is a model generator and solver 
interface used for linear, nonlinear and integer programming 
problems [Ref. 7]. GAMS generates mathematical models 
using algebraic statements which are easily read, easily 
modified and readily transported from one computer to another. 
The X-System [Ref. 6] is a solver used to obtain solutions to 
a wide range of mathematical programming problems including 
MIPS like the MCCES scheduling model. The version of the X- 
System used to solve this model resides on the Amdahl 5995- 
700A dual processor mainframe computer at the Naval 
Postgraduate School . 

A. PENALTY ASSIGNMENTS 

Critical to the behavior of an elastic optimization model 
is the proper assignment of penalty values for variables in 
the objective function. The penalties are assigned by the 
user and have values relative to the importance the user 
places on the violation of the constraints represented by the 
elastic variables. Additionally, the penalties must be 
proportional to the costs of other variables which appear in 
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the objective function. For instance, a large penalty 
associated with the elastic variable extracl^, along with a 
small value of the cost of a student in inventory, would 
result in placing a large number of students in inventory 
before a violation of CPTY^ (constraints (E5)) would occur in 
the model. Basic decision variable costs and elastic variable 
penalties used here in solving the MIP are contained in Table 
1 of the Appendix. 

B. SPECIAL CONSIDERATIONS 

Large integer programs are known to be computationally 
difficult to solve. For a fiscal year which has 236 training 
days and for a sample set of courses (17 different courses) 
from Communications Electronics Maintenance School (CEMS) , the 
elastic model generates a MIP which has 13,504 constraints, 
40,458 variables and 297,266 nonzero elements in the 
coefficient matrix. This is a relatively large integer 
program that takes 115 cpu seconds on the Amdahl to generate 
and could not be solved in a reasonable amount of time. 
Reducing the size of this model and capitalizing on special 
features of the X-System can lessen the difficulty of 
obtaining an optimal solution. 

1. Network Constraints 

An optional basis factorization in the X-System takes 
advantage of pure network constraints for improved efficiency 
[Ref. 8]. In the basic model formulation presented in 
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Chapter III, constraints (4) have a consecutive one's property 
for the variable z^i and a reformulation of these constraints 
into pure network constraints is accomplished by a standard 
procedure [Ref. 9] . The procedure requires the 

introduction of a slack variable, to convert constraints 

(4) from inequality to equality constraints. 

t 

^CT Vc,t^l (4') 

T=t- (INTVc-1) 



The following elementary row operations are then 
performed on constraints (4'): (a) the new first row, i = 1, 

is the same as the original first row, (b) for i = 2 ,..., last 
row, the new row i is obtained by subtracting original row i- 
1 from original row i. When the elastic variable, freqvio^i, 
is incorporated into the transformed constraints, constraints 
(E4a) and (E4b) result. 

t 

52 ^cT " fi^eqvio^t = 1 V c,t=l (E4a) 

T = t- (INTVc-1) 

^Ct ■ ^C(t-INTV_) '^ct “ 

(E4b) 

- freqviOj,^ + f reqvio^ (t-i) = ^ V c,tk2 

These new constraints are mathematically equivalent to 
constraints (4) but are pure network constraints, with each 
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variable having a +1 and a -1 coefficient in the combined 
constraints (E4a) and (E4b) . 



Constraints (5) also have the consecutive one's 
property for z^t- Since constraints (E4a) and (E4b) establish 
a maximal set of pure network constraints over the variable 
z^i, a reformulation of constraints (5) will not allow 
additional rows to be factored and is not likely to improve 
efficiency. 

2 . Planning Period Scaling 

The number of time periods in the planning horizon are 
identified as the number of training days in a fiscal year. 
The length of a class scheduled in the planning horizon is 
designated as the number of days necessary to complete a class 
of course c. All the decision variables and nearly all the 
constraints of the model are indexed by t, training days. To 
reduce the number of variables and number of constraints and 
thereby reduce the difficulty of obtaining an optimal solution 
to the MIP, all parameters associated with t are divided by 
the number of training days in a training week. 

For the MCCES scheduling problem, five training days 
are considered a training week and scaling of the sets indexed 
by t by the number of training days in a training week reduces 
the number of time periods in a fiscal year from 236 training 
days to 47 ( = f236/5l ) training weeks. References to all 
variables indexed by t are modified to reflect the change to 
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a weekly resolution. The length of a course is redefined as 
fLysl and represents the number of training weeks necessary 
to complete a class of course c. INTV^. is also scaled, 
flNTVc/sl , and redefined as the minimum number of training 
weeks between starting dates for course c. As a convention, 
reference to a training week will be taken to mean the first 
day of that training week. 

Planning period scaling has a dramatic effect on the 
elastic model by reducing the number of constraints and number 
of variables by nearly 80%, but some loss of resolution in the 
optimal solution occurs too. This is most easily seen by an 
example. Consider a simple sequence of two courses with no 
students in inventory and with the first course requiring 31 
training days to complete. In a model having daily 
resolution, a class of the first course which has 10 students 
and starts on training day one finishes on training day 31. 
In an optimal case, a class of the second course might start 
on training day 32 with 10 students in the class and require 
no students being placed in inventory. If the same two 
courses are scheduled by a model having weekly resolution of 
the time periods, a class of the first course with 10 students 
starting at the beginning of training week one finishes at the 
beginning of training week seven ( = [ 31 / 5 ] ) . The earliest 
possible time a class of the second course may start is the 
beginning of week eight since students can only start a class 
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of the second course on the next increment of time after 



completion of the first course. This schedule also does not 
require students to be placed in inventory in the model 
because placing students in inventory may only occur on the 
next increment of time after completion of a first course. 
Examination of the schedule obtained using a weekly time 
period resolution reveals a hidden accumulation of 40 man-days 
of inventory. Since the class of the first course finishes on 
the first day of training week seven, training day 31, the 10 
students must wait until the first day of training week eight, 
training day 36, which places each student in inventory for 
four training days: training days 32 through 35. The "slack" 
of a course c, slk^, is defined as the number of days before 
the end of a training week the course finishes and may be 
computed as 5 - MOD 5. The next section will discuss how 
slack and the resultant hidden inventory introduced by 
planning period scaling can, at least partially, be removed. 

The planning period scaling procedure reduces the size 
of the model generated from the data set identified earlier 
from 13,504 to 2,731 constraints, from 36,630 to 7,348 
variables and from 293,962 to 26,729 nonzero elements in the 
coefficient matrix. This is a significant reduction in the 
overall size of the model and the model generation time is 
reduced from 116 to eight cpu seconds. 
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A planning period scaling factor of five is a logical 
choice because a "normal" training week corresponds to the 
weekdays of a calendar week. Smaller scale factors could be 
used for greater accuracy but model size and solution times 
will increase. For an actual daily schedule (scale factor = 
1) , a solution could not be obtained using the elastic model 
after significant computational testing. 

C. SLACK REMOVAL AND CONVERSION TO A DAILY SCHEDULE 

Scaling of the time periods in the planning horizon by the 
number of training days in a week results in a schedule which 
specifies the starting dates of classes in terms of the week 
in which a class of a course will start. Using the convention 
adopted earlier, the conversion from week to day is simple: 
A class which starts in week t is synonymous with a class 
which starts on training day 5(t-l) + 1. 

For classes having starting dates rescaled to a daily 
resolution (t now denotes a training day) , the variable z^ 
equals 1 when a class of course c begins on training day t and 
equals 0 otherwise. (Remember that at most one class of a 
course can begin on a given day.) The hidden delay or slack 
caused by forcing classes to begin on the first day of a week 
is at least partially removed by the following procedure; 

• For each course c, compute slk^ from the weekly schedule 

• Convert the weekly schedule into a daily schedule 



35 



• For each course c , define echelon^ to be the echelon of 
course c where the echelon of the initial course is 
defined to be 1 {see Figure 1) 



• For all t = 1, . . . ,236 

For all courses c 

If a class of course c begins on day t 

Define the starting date of the class as stdtd= t 
If this is the first class of course c 
Define stdfcl^ = t 

• For e = 1 to number of echelons 

For each course c in echelon e 
For all courses c" € C® 

For all training days t = 1,...,236 
If 2,, = 1 

temps 1kg slkg 

While [ (stdtg,, > (stdfclg- + h^.) ) and 

(stdtg,, > (stdfclg. + INTVg-) ) and 
(stdtg., > 1 ) and 
(tempslkg >0) ] 

stdtg.j stdtg-i - 1 

tempslkg tempslkg - 1 
stdfclg. stdtg., 
slkg, slkg. + slkg 



Adjusting starting dates to remove the slack could cause 
constraint violations not already existing in the schedule. 
The logical comparison (stdtg.i > (stdfclg. + Lg.) ) ensures that 
CPTYg. is not exceeded by any amount greater than the weekly 
schedule. The logical comparison (stdtg-i > (stdfclg. + INTVgO ) 
ensures that no further violations of INTVg are committed and 
the logical comparison (stdtg-t > 1) ensures that the start date 
of a class remains in the fiscal year being scheduled. 
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D. NONINTEGER CLASS SIZES 



The fractional estimate of the attrition, ATTR,, is likely 
to cause noninteger class sizes in the schedule, so a 
modification to the flow balance constraints (constraints 
(E6) ) is made here. 



hc(t-D* Vc(t-L,)- NATTR^z^„.i.^,2 Vc't* Vc.tsl (e6>) 

c' € 



where NATTR^ is an integer estimate of the number of students 



who begin 


course 


c but 


fail 


to complete 


the course. 


Constraints 


(E6' ) 


require 


that 


the number 


of students 



scheduled to start a successor course on day t not exceed the 
number of students which have completed that course's 
predecessor. If MIN^ were equal to MAX^, a reasonable value 
for NATTRg would be [attr^MAX^ + . 5 J . However, two potential 
problems arise: (a) MAX^ > MIN^., which is always the case, and 
if MINj; students are scheduled for a class of course c, MIN^ - 
NATTRc might be a pessimistic estimate of those finishing (it 
is positive for this data, however), and (b) [atTR^MAX^ + . 5 J 
= 0 but ATTR^ 0. Case (a) is not a significant problem, at 
least for this data, because most classes scheduled are near 
MAXc in size. For case (b) , when ATTR^ is zero, constraints 
(E6') and (E6) are identical. When, due to rounding, NATTR, 
is zero and ATTR,. is not zero (MIN^ is always greater than zero 
by design) , no students fail to finish in the model when, in 
practice, there is some small attrition. A potential solution 
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to this problem is to increase the demand, TIP^, for course c 
by fTIP,ATTR,l , i.e., TIP^ TIP, + [tIPcATTR,! . This procedure 
was not attempted as it would require some judgement from an 
actual scheduler working with an actual schedule to evaluate. 
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V. COMPUTATIONAL RESULTS 



The solution to the scheduling problem faced by MCCES 
requires the determination of the starting date for each class 
of a set of courses along with specifying the number of 
students in each class. The MIP developed in the previous 
chapters is implemented to minimize the delay encountered by 
students while waiting for additional courses. To determine 
improvement made by using the MIP, a comparison is made 
between a class schedule developed by current scheduling 
methods (manual scheduling) and a class schedule developed by 
the MIP. 

A. MANUAL SCHEDULE 

To establish a reference to which comparisons may be made, 
an existing, manually created schedule is examined to 
determine the amount of delay encountered by students waiting 
for additional courses. Schedules developed by manual methods 
are not significantly different from year to year and the FY 
1993 schedule of classes is analyzed. To determine the number 
of man-days students spent waiting for courses of the 
Communication Electronic Maintenance School in FY 1993, a 
model instance is created for daily (scale factor = 1) 
scheduling of classes. Then, the class size variables, y^, 
and the binary variables indicating the start of a class, z^i. 
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are fixed to the values corresponding to the schedule 
published for FY 1993. Using the cost and penalty parameters 
contained in Table 1 and the course data of Table 2 , the 
elastic model identifies 7,819 man-days of delay in that 
schedule. This amount of delay sets a baseline to which any 
improvements to the schedule may be compared. 

The manually generated FY 1993 schedule of classes 
violates several constraints. The minimum class size of 
course RDRFC (MINr^rfc) is equal to six and two classes of the 
course RDRFC in the manually generated schedule, on training 
days 2 and 122, have five students scheduled. Additionally, 
the minimum interval which must be maintained between starting 
dates of classes of course c are violated several times. 
INTVrfc iNTVysRc violated three times and INTVgrrc 

INTVrdrfc both violated once. In light of these 

observations, the use of an elastic model which allows 
penalized violations of selected constraints is reassuring. 

B. AUTOMATED SCHEDULE 

To determine if the use of the elastic course scheduling 
model decreases the number of man-days students spend waiting 
for additional courses {improves the schedule) , the elastic 
model (constraints (E6') in place of constraints (E6) ) is used 
to generate a schedule of classes for the data of Table 2 but 
with the values of and INTVg scaled to weekly values, ie. 
divided by five and rounded up if fractional. The cost and 
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penalty parameters of Table 1 are also used as inputs to the 
model . 

After the elastic model determines a schedule, the 
starting dates of courses are rescaled to a daily resolution 
and are then adjusted by the procedure presented in Chapter IV 
to reduce the slack. Then, the class scheduling variables, 
Zct/ and class size variables, y^t, are held constant and 
another instance of the elastic model for daily scheduling 
(scale factor = 1) is created and solved to determine the man- 
days of delay present in the schedule with slack removed. The 
solution to this instance of the model is obtained in a total 
of 374 cpu seconds. This time reflects generating and solving 
a weekly model (220 seconds) , slack removal and conversion of 
start dates (insignificant) , generating and solving of a daily 
version of the model (154 seconds) . The schedule specified by 
these steps has 2,964 man-days delay which is a decrease of 
4,855 man-days of delay from a manually developed schedule. 
The schedule developed by the elastic model meets all 
constraints except for a violation of the minimum class size 
for the course GRDR on training days 95 and 135. The minimum 
class size for GRDR is six and the number of students 
scheduled for the classes on training days 95 and 135 is five. 

To examine the behavior of the MIP to various input 
parameter values, several experiments have been performed. 
The parameters used in these experiments are identified in 
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Table 3 and the results of the experiments are summarized in 
Table 4. These experiments indicate that student delay is 
sensitive to penalty values but that solution times are only 
slightly affected. These experiments also indicate that some 
"tuning" of penalties may be necessary in practice, but that 
such tuning is not likely to cause computational problems. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



Comparison of the student delay in schedules constructed 
by manual and automated scheduling methods shows that a mixed 
integer programming approach can provide an improved schedule. 
For the set of courses examined, the amount of time students 
spend waiting for classes of a subsequent course can be 
reduced by about 4,855 man-days or over 62%. (This is about 5% 
of the total training man-days.) For the group of courses 
tested, the average delay per student is 11.7 training days 
(man-days delay/number students entering the initial course) 
under the manual schedule. For the automated schedule, the 
average delay per student is reduced to 4.4 training days 
(actual delay may be higher if NATTR, is significantly 
different from the actual attrition and if many classes are at 
or near their minimum class size) . This means that if an 
automated scheduling procedure is used, students attending 
Communications -Electronic courses could proceed to their 
subsequent assignments much more quickly. 

An alternate way to quantify the savings provided by an 
integer programming approach to scheduling is by assigning a 
dollar value to the delay. If the average daily wage of a 
student who is waiting for an additional course is $30 
[Ref. 10] , the schedule produced by the MIP would 
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save approximately $145,650 in salary of students waiting for 
subsequent courses. 

When the students in inventory are collected among all 
course groups at MCCES, the population may at times be quite 
large. Another saving created by reducing the number of 
students in inventory is the reduction in the number of 
supervisory personnel required to manage these students. 

The use of a mixed integer programming approach to the 
scheduling of classes at MCCES holds the potential for a large 
improvement in efficiency of the training of students at 
MCCES. The user should be aware that NATTR^ is typically a 
pessimistic attrition number based on class sizes being near 
MAX^. If class sizes scheduled are small, it might be 
reasonable to adjust NATTR^ downward. Reformulations of the 
models presented in this thesis using fractional attrition, 
ATTR,, may produce solutions with integer class sizes and such 
reformulations could be explored. Additionally, the 
development of a user-friendly interface would be valuable. 

The use of this type of scheduling model may be applicable 
with little or no modification to many armed forces schools. 
This model should apply to other schools which have multiple 
instances of sequential courses and which train large numbers 
of students each year. 
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APPENDIX 



MODEL INPUT PARAMETERS 
TABLE 1. COST AND PENALTY VALUES 



Parameter 


Value 


K 


0.01 




-0.05 


MAXPEN^ 


8.0 


MINPEN^ 


0.8 


TIPPEN^ 


8.0 


tippen; 


8.0 


FREQPEN, 


80.0 


EXTRAPEN, 


80.0 



The values of parameters in this table are the values 
which are used as inputs to the elastic model presented in 
Chapter III. These values may be changed by the user to 
achieve different behavior from the elastic model. 
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TABLE 2. FY 1993 DATA FOR COMM-ELEC MAINTENANCE SCHOOL 



COURSE 


MAX^ 


MIN^ 


Lc 


TIP, 


INTV, 


CPTY, 


ATTR, 


RFC 


30 


15 


98 


668 


5 


12 


.169 


RDRFC 


12 


9 


20 


49 


20 


1 


.030 


MRC 


12 


1 


98 


9 


12 


1 


.037 


TSRC 


14 


8 


71 


98 


30 


4 


.075 


TACC 


12 


5 


98 


15 


98 


2 


.086 


TDCC 


9 


6 


98 


5 


98 


8 


0 


ACEOC 


8 


1 


98 


12 


69 


1 


0 


RFC 


24 


12 


30 


359 


9 


4 


.015 


AIMS 


24 


8 


50 


38 


50 


1 


.067 


GRDR 


8 


8 


37 


12 


37 


1 


.317 


TAOM 


6 


3 


l<i 


12 




1 


.027 


ARRC 


12 


4 


98 


23 


17 


8 


.086 


MEMC 


15 


12 


100 


30 


50 


8 


0 


GRRC 


24 


12 


78 


292 




6 


.072 


ARDRA 


12 


4 


98 


14 


69 


1 


0 


ARDRB 


12 


4 


40 


17 


40 


1 


0 


PSMC 


8 


4 


17 


18 


17 


1 


0 



The data in this table is specific to the set of courses 
for the Communication- Electronics Maintenance School. The 
values are required as inputs to the elastic model of Chapter 
III. 
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TABLE 3. EXPERIMENTAL COST AND PENALTY VALUES 



Parameter 


Experiment 


Experiment 


Experiment 


Experiment 




1 


2 


3 


4 


K 


0.01 


0.01 


0.01 


0.01 


Be 


-0.05 


-0.05 


o 

o 


-0.05 


MAXPEN„ 


8.0 


4.0 


2.0 


1.0 


MINPEN^ 


0.4 


0.4 


0.2 


0.1 


TIPPEN^ 


8.0 


4.0 


2.0 


1.0 


tippen; 


8.0 


4.0 


2.0 


1.0 


FREQPEN, 


80 . 0 


40.0 


20.0 


10.0 


EXTRAPEN, 


80.0 


40.0 


20.0 


10.0 



The values of the parameters of this table are variations 
of the parameters contained in Table 1. The values are 
changed to examine the response of the elastic model to 
different input parameter values. 
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TABLE 4. EXPERIMENTAL RESULTS 



Experiment 


Solution 
Time (sec) 


Man - days 
Delay 


Violations 


1 


374 


2964 


yORDR.95 = 5 
yCRDR.135 = 5 


2 


387 


3416 


yCRDR.135 = 4 


3 


410 


4688 


yoRDR.eo = 4 


4 


429 


3248 


yCRDR.70 = 3.75 

yARDRB.51 = 3 



This table displays the response of the elastic model to 
the differing parameter values of Table 3. In addition to the 
violations identified for Experiment 4, five other classes in 
the schedule have fractionated class sizes. 
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